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ABSTRACT 



We have compiled a sample of 3041 spiral galaxies with multi-band gri imag- 
ing from the Sloan Digital Sky Survey (SDSS) Data Release 7 and available 
galaxy rotational velocities, V, derived from Hi line widths. We compare the 
data products provided through the SDSS imaging pipeline with our own pho- 
tometry of the SDSS images, and use the velocities, V, as an independent metric 
to determine ideal galaxy sizes (R) and luminosities (L). Our radial and lumi- 
nosity parameters improve upon the SDSS DR7 Petrosian radii and luminosities 
through the use of isophotal fits to the galaxy images. This improvement is 
gauged via VL and RV relations whose respective scatters are reduced by ~8% 
and ~30% compared to similar relations built with SDSS parameters. The tight- 
est VRL relations are obtained with the i-band radius, i?23.5,i, measured at 23.5 
mag arcsec"^, and the luminosity I/23.5,i, measured within i?23.5,i- Our VRL scal- 
ing relations compare well, both in scatter and slope, with similar studies (such 
comparisons however depend sensitively on the nature and size of the compared 
samples). The typical slopes, b, and observed scatters, a, of the i-band VL, RL 
and RV relations are bvL = 0.27 ± 0.01, buL = 0.41 ± 0.01, bRv ^ 1.52 ± 0.07, 
and avL = 0.074, grl = 0.071, any = 0.154. Similar results for the SDSS g and 
r bands are also provided. Smaller scatters may be achieved for more pruned 
samples. We also compute scaling relations in terms of the baryonic mass (stars 
-I- gas), Mbar, ranging from M^ar — 10^'^ Mq to 10^^-^Mq . Our baryonic velocity- 
mass (VM) relation has slope 0.29 ± 0.01 and a measured scatter ameas — 0.076 
dex. While the observed VL and VM relations have comparable scatter, the stel- 
lar and baryonic VM relations may be intrinsically tighter, and thus potentially 
more fundamental, than other VL relations of spiral galaxies. 
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Subject headings: galaxies: dynamics — galaxies: formation — galaxies: kinematics 

— galaxies: spirals — galaxies: structure — dark matter 
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1. Introduction 

Fundamental scaling relations for spiral galaxies are known to emerge from the 
combination of observed galaxy rotation velocity, V, total luminosity L, and size, R. For 
instance, the VL relation, also known as the TuUy-Fisher relation (TuUy & Fisher 1997), 
probably defines the fundamental plane of spiral galaxies. That is, the scatter of the VL 
relation cannot be reduced by considering any other third parameter (e.g., Courteau & Rix 
1999, hereafter CR99; Courteau ct al. 2007, hereafter C07; Button et al. 2007, hereafter 
D07; Pizagno et al. 2007). The study of galaxy scahng relations also enables a direct 
comparison with theoretical models of galaxy formation (e.g., Pizagno et al. 2005; D07; 
Avila-Reese et al. 2008, hereafter AR08; Dutton & van den Bosch 2009; Button et al. 2011, 
hereafter Bll). These and other studies suggest that the simultaneous matching of the VL 
and RL relations, whilst matching the observed galaxy luminosity function and reproducing 
the shape of galaxy surface brightness profiles is a challenging task. For example, for 
standard disk galaxy models (e.g., Mo, Mao & White 1998; B07) with standard stellar 
initial mass functions (i.e., Kroupa 2001/Chabrier 2003) to match basic galaxy scaling 
relations, B07 and Bll showed that halo expansion is required. This may be reahzed 
through dynamical friction on baryonic clumps and/or supernova driven mass outflows 
(e.g., Navarro, Eke, & Frenk 1996; El-Zant, Shlosman, & Hoffman 2001; Mo & Mao 2004; 
Governato et al. 2010; Cole, Behnen & Wilkinson 2011). 

It has also been stated that the luminosity of a spiral galaxy is a poorer tracer of its 
circular velocity than baryonic mass, -Mbar (McGaugh et al. 2000; McGaugh 2005). The 
latter is defined as the sum of the luminous mass, M*, and the gas mass, Mgag. is 
usually obtained by multiplying the total extrapolated luminosity, measured in a specific 
wave band, by a suitable stellar mass-to-light {M/L) ratio. Mgas is measured directly from 
the Hi fiux, using a correction factor of 1.4 to account for the mass fraction of helium. At 
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low total galaxy mass, Mgas can be a significant fraction of M^,, thus raising the question 
whether L, M^, or Mbar is a better match to V. The latter, which is also known as the 
"baryonic Tully-Fisher" or "BTF", may also have a different slope than the standard VL 
relation (Bell & de Jong 2001; Verheijen 2001; McGaugh 2005; Gurovich et al. 2010). The 
samples that have been used for BTF studies have typically included fewer than ~ 50 
galaxies based largely on multiple, heterogeneous samples. For their sample of 243 galaxies, 
McGaugh et al. (2000) found that the BTF relation is more "fundamental" than the VL 
relation; deviations from the BTF relation may however exist (McGaugh & Wolf 2010). 
However, a BTF relation based on a large, statistical sample of galaxies is still lacking. The 
discrepancies between published BTF slopes and scatters also motivate a new study with 
as large a galaxy sample with accurate rotation velocities and Hi fluxes as possible. 

Interest in galaxy scaling relations also stems from wanting accurate distance 
estimators, which in turn is obtained via the suitable pairing of a distance-dependent and 
distance-independent galaxy parameters such as size, luminosity or colour with circular 
velocity. The infrared VL relation has typical distance errors of 15 — 20% (e.g., Aaronson 
et al. 1979; Pierce & Tully 1988; Gavazzi et al. 1999; Masters et al. 2006; hereafter 
M06; Saintonge & Spekkens 2011; hereafter SSll). Very large data samples and carefully 
measured galaxy parameters can reduce sampling error in these studies. The Sloan Digital 
Sky Survey (Abazajian et al. 2009; hereafter SDSS) is currently the largest data base of 
galaxy structural parameter^ It is thus relevant to ask if the SDSS library of galaxy 
scaling parameters yields the tightest possible scaling relations. 

A main goal of this paper is indeed to investigate the quality of galaxy scaling relations 
based on size, luminosity, and colour derived from SDSS data products. The intent is to 

-^The NYU Value-Added Galaxy Catalogue by Blanton et al. (2005) also offers a cross- 
matched collection of galaxy catalogs, using the SDSS library as a core. 
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compare SDSS pipeline data products with similar measurements extracted from SDSS 
images but using independent data reduction methods. 

We develop our analysis of VRL relations through two specific channels: i) we first 
compare the multi-band data products from the SDSS DR7 with our own galaxy structural 
parameters extracted from SDSS FITS galaxy images, and ii) we determine which of the 
VRL parameters yield the tightest luminous and baryonic scaling relations. Limitations 
of the SDSS data products will be addressed along the way (see also Masjedi et al. 2006; 
Lauer et al. 2007; Fathi et al. 2010). 

To set the foundations of the VRL relation, we follow Courteau et al. (2007, hereafter 
C07) whose sample comprised 1300 late-type galaxies. The sample that we present here, 
having more than ~ 3000 galaxies, is a two-fold increase over C07. We follow most, though 
not all, of the reduction methods and parameter corrections from C07. For instance, 
C07 used photometry and rotational velocities from four separate sources (Mathewson 
et al. 1992; TuUy et al. 1996; Dale et al. 1999; Courteau et al. 2000). A significant 
improvement of this study over C07 is our use of strictly homogeneous data for both the 
multi-band photometry (SDSS) and fine widths (S05/S07). C07 also used disk scale lengths 
as a measure of galaxy size in order to facilitate comparisons with galaxy formation models 
(D07). We consider disk scale lengths here as well, but will also extract other radial metrics 
that yield tighter VRL relations (see also SSll). 

The available SDSS parameters that are of interest to us, namely the galaxy size, 
R, and luminosity, L, are both distance-dependent. Thus, in order to test which of our 
measurements or the SDSS parameters yield the tightest scaling relations, we must compare 
each value against an independent metric which we take here as the distance-independent 
galaxy velocity. 

Our study benefits from the availability of rotational velocities, V, from spatially- 
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unresolved neutral hydrogen (Hi) 21 cm spectral line widths, for nearly 9000 spiral galaxies 
(Springob et al. 2005; hereafter SOS; Springob et al. 2007, hereafter S07). Cross-correlation 
of the SOS and S07 line width catalogs with the SDSS will define our target sample. 

Our paper is organised as follows: In Section [2| we cross-correlate the SOS and S07 Hi 
line width catalogs with the SDSS DR7; this yields a sample of 3041 disk galaxies with 
accurate photometry and line widths. We discuss the extraction of radial and luminous 
parameters from SDSS images in Section |3j Final galaxy parameter corrections are applied 
in Section |4] and the data table of galaxy structural parameters is presented in Section |5} 
In Section [6} we compare the SDSS and our galaxy parameters against an independent foil, 
here chosen to be the (distance-independent) rotational velocity, V, from S0S/S07. We 
determine the "best" scaling parameters with which to build the tightest VRL relation in 
Section [7] and present, in Section |8} the largest baryonic Tully-Fisher sample to date. A 
summary of our results is presented in Section [9] 



2. Data Sample 

In order to simultaneously test the reliability of SDSS data products and establish the 
most comprehensive VRL scaling relations to date, we seek not only the largest but the 
most homogeneous compilation of galaxy structural parameters to date. This can be done 
through the cross-correlation of the large compilation of ~9000 galaxies within cz < 28, 000 
kms~^with Hi line widths by SOS and S07 with the multi-band ugriz photometry provided 
by the SDSS. 

The SDSS archive query with S0S/S07 targets yielded FITS images for 4260 galaxies. 
That sample was examined visually to eliminate galaxy pairs, interacting galaxies and 
images with bright foreground stars, and any galaxy whose light profile could not be 
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extracted (e.g., very faint targets). Due to the survey nature of the SDSS, target galaxies 
may often appear close to the edge of the CCD image frame and subsequently compromise 
the photometric analysis. A large fraction (~ 40%) of the more luminous galaxies in our 
sample which suffer from these "edge effects" had to be discarded. 

We were left with a final sample of 3041 galaxies for which both S05/S07 rotational 
velocities and acceptable SDSS FITS images are available. We have subdivided this large 
sample, dubbed "Sample A", into three subsets (Samples B, C, D) to investigate the effects 
of distance uncertainties and inclination on galaxy scaling parameters. 

The main sample and subsets include the following systems: 

(A) Full Sample - 3041 galaxies from S05/S07 found in SDSS DR7 with imaging suitable 
for surface photometry; 

(B) Best Inclinations - 1725 galaxies from Sample A with moderate inclinations 

40° < i < 75°. This removes the nearly face-on galaxies that would suffer from large 
uncertainties in their deprojected rotational velocity, as well as the more inclined 
galaxies (i > 75°) whose disk may be significantly obscured by line-of-sight dust and 
whose surface area for isophotal photometry is barely visible; 

(C) Best Distances - 1076 galaxies in Sample A with the best distance determinations from 
high quahty spectroscopic data of S07 (referred to as "SFI-I— 1-" by S07) ; 

(D) Best Inclinations and Distances - 652 galaxies from the intersection of Samples B and 
C. Sample D is thus our "Best Sample". 
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3. Light Profile Extraction 

3.1. SDSS FITS Image Photometry 

We extract light profiles from the SDSS gri images of all the Sample A galaxies using 
the surface brightness profile extraction methods presented in Courteau (1996) and adapted 
to SDSS images by McDonald et al. (2009; 2011). Given the lower signal-to-noise of the 
SDSS u and z bands (Blanton et al. 2001), our analysis will rely strictly on the gri bands. 

The image processing software XVISTJij^is the backbone of our surface brightness profile 
calculations. The XVISTA command PROFILE was used to fit azimuthally-medianed elliptical 
isophotes through the galaxy 2D light distribution. The position angle and ellipticity of 
each isophote could vary while the centre was kept fixed. The ellipticity which best defines 
the stable outer disk was visually determined and the isophotal solution was extended to 
faint light levels using that ellipticity to extract the deepest profile possible. The surface 
brightness was calculated from the median sky-subtracted fiux level over an elliptical 
contour. Intervening stars were masked and isophotal twists due to bright spiral arms were 
also smoothed out by hand. 

The higher signal-to-noise of the z-band SDSS images, and the lesser sensitivity of 
galactic dust at redder wavelengths, make the i-band the filter of choice for galaxian 
surface photometry. The profile templates (position angle and ellipticity) for each galaxy 
is based on the i-band image and then applied to both g and r images. This ensures that 
galaxy colour profiles are extracted over the same physical regions of a galaxy. The surface 
brightness levels at each pixel were computed in each band and scaled using the photometric 
zero-point, airmass and extinction coefficients supplied in the SDSS-drObj photometric 
calibration files. The surface brightness profiles for the 3041 galaxies in Sample A were all 



^see http : / / ganymede . nmsu . edu/holtz/ xvista/ 



- 10 - 



inspected individually for quality control. 

3.2. SDSS Sky Measurements 

The galaxy light profiles were all sky-subtracted; this operation being clearly the 
largest source of error in surface brightness extraction. The SDSS image headers already 
include an estimate of the sky background level. However, the sky backgrounds for bright 
extended galaxies in the SDSS DRl - DR7 are known to be systematically over-estimated 
resulting in the under-estimation of galaxy fluxes and sizes by as much as 20-30% (Masjedi 
et al. 2006; Lauer et al. 2007; Abazajian et al. 2009; McDonald et al. 2011). Conversely, 
for fainter systems such as LSB galaxies. West et al. (2010) report an over-subtraction of 
the sky level in the SDSS photometric pipeline. We now present our investigation of sky 
subtraction uncertainties based on DR7 SDSS z-band galaxy image^ 

The SDSS sky background level, skysDss, is measured from an initial estimate of the 
median pixel value of every fourth pixel in the image clipped at ~2.3 a, where a is the rms 
noise of the image. The brightest stars, but not necessarily their wings, are thus removed 
from the image and the sky level is recomputed to provide a final background estimate. 
The different sky levels as provided by SDSS for Sample A galaxies are shown in Figure [T] 
for the three SDSS gri bands. 

To investigate the deviations in sky measurements, we compute two other sky 
measurements from the SDSS i-band FITS images: i) skyFuu, is the median pixel value over 
the full image frame and, ii) skyxrue is the median sky background level from five boxes 

^Note that background subtraction improvements in DR8 have resulted in more reliable 
photometry of large galaxies (Blanton et al. 2011). The treatment about sky uncertainties 
presented in our paper is still relevant for smaller galaxies, whether in DR7 or DR8. 
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placed interactively away from the galaxy image. The latter was computed for a sample of 
30 galaxies with the largest and smallest angular sizes. In principle, skyxrue provides the 
most realistic estimate of the sky background. 

We measure the difference Asky in i-band sky levels with respect to skyxme as, 



skyi-skyTrue 

Askyi = . 1 

skyTrue 

Figure |2] shows the results for Askysoss and Askypuu for our subsample of 30 selected 
galaxies. First, we note that the largest deviations in both panels are for galaxies whose 
angular diameter exceeds 7'5. Fortunately, most of our galaxies have much smaller sizes. 

The top panel of Figure |2] confirms that the SDSS sky levels are overestimated with 
respect to our "best" estimates, sky True- The percent difference for most small galaxies is 
however less than ~0.2%. The offset between the "Full Sky" estimate and our best sky 
measurement (bottom panel) is surprisingly comparable to that with the SDSS sky level 
(top panel). In most cases, the deviations are less than 0.4%. 

The effect of over- and under-subtracted skies on a typical light profile is shown in 
Figure |3] for sky background errors at the 0.2% (light blue points), 0.5% (dashed line) 
and 1.0% (dark solid line) levels on the i-band surface brightness profile of UGC 5651 (an 
intermediate-size Sc galaxy centered on the image frame). The dashed black line marks the 
H = 23.5 mag arcsec"^ level and it can be seen that most isophotal measurements above 
the /i = 23.5 mag arcsec"^ level are largely free of normal sky subtraction errors. We will 



compute in Section 6.2 the effect of sky subtraction errors on measured sizes and fluxes. 



Based on these tests, we have adopted the SDSS sky levels, skysDss, for the sky 
background subtraction of our surface brightness profiles. Caution must still be taken when 
utilizing SDSS sky background estimates, especially for large galaxies. 
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Figure |4] shows the median surface brightness errors versus median surface brightness 
in the gri bands for all the galaxies in Sample A. Our sky-subtracted galaxy light profiles 
are quite stable and reliable down to surface brightnesses of ;U ~ 27 mag arcsec"^, where 
surface brightness errors typically do not exceed 0.4 mag arcsec"^. 



3.3. Structural Parameter Extraction 

We now present the measurements of structural parameters in the gri bands for the 
3041 galaxies in Sample A. We seek isophotal, effective and physical size and luminosity 
parameters for the construction of VRL relations. 

Isophotal parameters are measured at a specific surface brightness level, effective 
parameters are measured at a radius which encloses half of the total light, and physical 
sizes can refer to a disk scale length or an exact radius in kpc. 

The total luminosity of a galaxy can be estimated by extrapolating its light profile 
to infinity assuming an extended exponential disk. Given that the SDSS photometry is 
already fairly deep (Figure |4]), the profile extrapolation adds only a few percent of the light 
to the total. The exponential disk model is given by 

I{R) = Ioexp{-R/Ra), (2) 

where /q is the extrapolated central surface brightness in counts and R^ is the scale length 
of the exponential disk. In magnitude units. Equation ^ becomes 

^(i?) =/io + 1.086(i?/i?d), (3) 

where the units of fi are mag arcsec"^. The disk light, L^, in counts, is the cumulative 
surface brightness under the exponential fit from r = to infinity, 

Ld= / I{r)dr. (4) 

^0 
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Disk extrapolation fits are highly susceptible to the fit baseline and the presence of 
spiral arms. Indeed, disk scale lengths {Rd) and central surface brightnesses (yUo) may vary 
by as much as 20% from author to author (Knapen & van der Kruit 1991; Courteau 1996; 
Courteau et al. 2011). Here we adopt the interactive "marking the disk" technique to fit 
an exponential function to the disk surface brightness profile (Giovanelli et al. 1999; COT). 
The inner and outer radii for the fit baseline encompass the region of the SB profile that is 
dominated by the disk. Full 2D bulge-disk model decompositions of the SDSS images will 
be reported elsewhere; for our current purposes, the "marking the disk" technique is fully 
adequate. 

Besides the measurement of a disk scale length and disk central surface brightness, we 
can also compute from these fits the isophotal radius, -R23.5, and apparent magnitude, m23.5, 
within the 23.5 mag arcsec"^ isophote. This straightforward measurement does not depend 
on any disk fitting or parameterisation. Other isophotal radii can be measured but -R23.5 



can be shown to yield smallest scatter of the RL relation (e.g., Courteau 1996; Section 6.2) 



We also compute the half-light (or effective) radius Re containing 50% of the total 
extrapolated light. The effective surface brightness is defined as fie = fJ'{Re)- We also define 
other parameters: mext is the total magnitude of the extrapolated light profile and m2,2 the 
magnitude integrated within the radius R = 2.2i?d- The latter corresponds to the peak 
of the rotation curve of an idealized pure exponential disk (e.g.. Freeman 1970; Binney & 
Tremaine 1987). 

We now have various sets of radii and apparent magnitudes based upon our extrapolated 
light profiles for 3041 galaxies in g, r and i-bands. Colours can also be computed from 
various magnitude combinations. In Section |6] we quantify the robustness of these structural 
parameters in the context of the VRL scaling relations. 

Our analysis benefits from the independent reduction of 211 galaxy light profiles by two 
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of us (MH and YZ). The light profiles for nine of these galaxies are presented in Figure |5] 
as blue (reduced by MH) and red (reduced by YZ) profiles. Comparing these profiles shows 
that differences between MH and YZ, if any, are small and purely random, despite slightly 
different treatments in the few interactive aspects of the SB profile extraction. These two 
treatments differ only in the smoothing of isophotal twists and luminosity extrapolation 



of the outer disk (see Section 3.2). MH extracted the light profiles for all 3041 galaxies in 
Sample A and we adopt her final catalog. The extensive comparisons with YZ reinforce the 
reliability of our sample. 



4. Final Data Corrections 

We now correct the extracted structural parameters (e.g., sizes, luminosities, surface 
brightnesses) following the prescription of C07 but now accounting for the SDSS filter 
system. The velocities derived from the Hi line profiles are corrected according to the 
prescriptions of SOS and S07. The inclination of each galaxy is measured from the ellipticity 
of the outer disk isophotes of the SDSS i-band image and corrected for projection effects 
using Holmberg's oblate spheroid description, 

where the ellipticity, e = 1 — 6/a, go = c/a is the axial ratio of a galaxy viewed edge-on, a 
and h are the disk major and minor axes, and c is the polar axis. 

In order to determine the variation of projected ellipticity with colour or morphological 
type, we have used the Galaxy Zoo classification (Lintott et al. 2008). We selected all of 
our galaxies classified as edge-on by > 80% of all the Galaxy Zoo survey participants (and 
eliminated all those classified as ellipticals or merger by > 90% of them) yielding a sample 
of 1377 candidate edge-on galaxies. We then visually inspected the SDSS images for all 
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these galaxies and eliminated a further 506 systems deemed not perfectly edge-on. For 
each of the remaining 871 edge-on galaxies, an ellipticity was computed using the adaptive 
second moments, e+ and tx-, provided by the SDSS pipeline. From these moments, the 
adaptive ellipticity was computed following Kautsch (2009). Finally, this ellipticity was 
converted to an axial ratio (go = 1 — e). 

Figure |6] shows the variation of SDSS-derived axial ratios in the i-band versus mean 
galaxy colour, computed at the effective radius, for 871 edge-on galaxies each depicted 
as a black dot. The morphological correspondance (top axis) was derived from SDSS 
colours tabulated in McDonald et al. (2011). The red dots are median averages per 
given morphological bin. The trend of go naturally increases with redder, progressively 
bulge-dominated galaxies. The go distribution is relatively flat for blue (late-type) galaxies, 
with an upturn at (7 — z ~ 1.3, corresponding to Sa/SOa galaxies. Several SO galaxies have 
go ^ 0.8, reminiscent of the Sombrero galaxy. The few late-type systems with go > 0.5 
are likely a failure in the adaptive moment pipeline. Note that similar go distributions are 
obtained for the g and r bands. 

Based on Figure [6} we adopt go — 0.13 as the minimum axial ratio, and thus 
intrinsic thickness, of spiral galaxies. This value agrees with a similar study by Giovanelli 
et al. (1994), though a value as high as 0.2, as reported by Lambas et al. (1992), is also 
realistic considering the intrinsic errors on galaxy radii. For their study of the baryonic 
TuUy-Fisher relation (see Section [s]). Stark et al. (2009) adopted go = 0.15. The effect of 
using go = 0.13 or 0.2 is inconsequential for our study. 

The apparent ((yri-band) magnitudes have been homogeneously corrected for internal 
extinction Aj, external galactic extinction Ag and /c-correction according to 

m = mobs - Ai-Ag- A^. (6) 
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The inclination-dependent internal extinction correction in magnitudes is given by 

A = 7log(V«) (7) 

where 7 is interpolated in the g, r and i bands from the corrections of TuUy et al. (1998) 
given the corrected line widths logl^ from SOS, 

73 = 1.51 + 2.46(logM^-2.5), (8) 
7r = 1.25 + 2.04(logVF- 2.5), (9) 
7i = 1.00 + 1.71 (log ly- 2.5). (10) 

The Galactic extinction, Ag, from Schlegel et al. (1998) is provided by the SDSS. The 
cosmological /c-correction, A^, is calculated from the template of Blanton & Roweis (2007) 
for each SDSS filter. 

The absolute gri magnitudes, M, are calculated in their respective band as 

M = m - 5 log + 25; Dl = VcmbH^^ (11) 

where the Hubble constant Hq = 71 km sec"^ Mpc"^ (Komatsu et al. 2011) and is 
the distance in kpc derived from the velocities in the rest frame of the cosmic microwave 
background, Vcmb, as compiled by S05 and S07. For cahbration to solar units, we use 
-^©,3=5.12, M0^r=4.68 and M0_j=4.57 (York et al. 2000); we use in the conversion from 
apparent to absolute quantities. 

We correct both effective (/ie) and central (/^o) surface brightnesses for internal and 
external extinction following C07: 

III ^ lie + 0.5 log(a/6) -Ag- 2.5 log(l + zf, (12) 
n'o^Ho + 0.5 log(a/6) - Ag - 2.5 log(l + zf. (13) 
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These are reported in a data table (Table [9] as presented Section [5]) but not used for further 
analysis in this paper, except in Figure [221 Radii measurements may be corrected following 



Giovanelli et al. (1995), as also reported in C07. However, we shall find in Section 7.2 that 
such corrections do not yield a reduction of the RL and RV relations scatter. In view of the 
uncertain nature of most radial scale corrections, our analysis relies on uncorrected radii. 



4.1. Parameter Uncertainties 

We follow M06, C07, AR08, SSll and many other studies of galaxy scaling relations 
in estimating the uncertainties for each scaling relation parameters logV^, logi? and logL. 
Uncertainties in logL assume 10-15% errors for each of the parameters 6/a, 7, Agai, and 
D in Equation (|6]). That is, 

^logL - ^mag + + +f^log(b/a) + ^Ag^i + ^Afc + ^Z)> (14) 

where cxmag is the statistical error in the raw magnitude measurement. A more complete 
assessment of VRL parameter errors is found in SSll. Our results are however practically 
the same. 

Typical measurement errors on L amount to 20%, or dlogL = 0.09. This value is also 
obtained by M06 and SSll for their measurements of the absolute luminosity, M/, in the 
J-band. The uncertainties on the rotational velocities as compiled by S05/S07 amount to 
~5%, or dlogV = 0.02 (SSll use dlogV^ = 0.03). We also assume a 15% uncertainty, or 
dlogRd = 0.07 on disk scale lengths (e.g., MacArthur et al. 2003; Courteau et al. 2011) and 
a 7% error, or dlog-R=0.03, on isophotal radii (same as SSll). 

In the VRL fits that follow (Section [T]), the quantities dlogL and dlog\^ are computed 
independently for each galaxy while dlogi? for various radii measurements (-R23.5, Rd, Re) 
is set to the mean values above for each galaxy. Using this method or assigning a constant 
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error to each parameter plays no role in the final scaling relations. 



5. Data Table 

We now present our extensive compilation of galaxy structural parameters. The 
recessional and rotational velocities, Hi flux, and morphological classiflcation are from 
S05/S07. All other structural parameters have been derived by us from SDSS gri images 



as described in Section 3^ The distribution of several of these parameters is shown in 
Figure [7] for the full Sample A in red and for our Sample D with best constrained distances 
and inclinations in blue. Figure [7] shows a broad parameter coverage, though none of which 
can be deemed complete. The parameter distributions are roughly (log-)normal with an 
emphasis on Sb/Sc morphological types. 

Table [9] gives a list of structural parameters for all the Sample A galaxies. The first 
few entries of the catalog are shown here; please see the electronic version for the full list. 
The entries are arranged as follows: 



Col. (1): UGC or AGC galaxy name; 

Col. (2): i, the corrected galaxy inclination in degrees using = 0.13 in Equation ([s]); 

Col. (3): VcMB, the recession velocity of the galaxy relative to the cosmic microwave 

background frame in kms~^ as calculated by S05 and SOT. The S07 estimates of Vcmb 
supersede SOS in the case of multiple measurements; 

Col. (4): Vj-ot, the corrected rotational velocity in kms~^ from the Hi line width profiles 
measured at 50% of total fiux by SOS & S07 (PVpso)- S05/S07 line widths are corrected 
for instrumental effects and redshift broadening. V^ot is corrected for the inclination 
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of the disk and distance but not for turbulent motion in the Hi gas. Sample D uses 
only rotational velocities from S07; 

Col. (5): T, the RC3 morphological type as listed by SOS & S07; 

Col. (6): Re, the half-light radius based on the total extrapolated z-band luminosity of the 
galaxy in arc seconds; 

Col. (7): R23.5, the radius at the z-band 23.5 mag arcsec"^ isophote in arc seconds; 

Col. (8): Rd, the disk scale length of the galaxy in arc seconds; 

Col. (9): m^^^, the total apparent magnitude of the galaxy including the extra light from 
the disk extrapolation to infinity; 

Col. (10): m23 5, the apparent magnitude within the i-band 23.5 mag arcsec"^ isophote; 

Col. (11): ni2 2, the apparent magnitude within 2.2i?d of the surface brightness profile. 
Non-entries in this column are due to a failed disk extrapolation, described in 



Section 6.3 



Col. (12): Hq, the corrected central surface brightness in mag arcsec"^; 

Col. (13): yUg, the corrected effective surface brightness measured in mag arcsec"^; 

Col. (14): g — r, the colour of the galaxy calculated from the difference in g and r 
extrapolated magnitudes rriext] 

Col. (15): g ~ i, the colour of the galaxy from the g and i bands as above; 

Col. (16): C28,j, the galaxy light concentration, C28 = 5 \og{Rso / R20) , where R20 and Rso 
enclose 20% and 80% of the extrapolated i-band light, respectively; 



- 20 - 



Col. (17): log Mhi, the neutral hydrogen gas mass in Mq derived by SOS from the 
corrected Hi hne flux measurements; 

Col. (18): logM^,, the stellar mass of the galaxy in units of Mq from the extrapolated 

i-band luminosity L'^^^^ ; times the stellar mass-to-light ratio M^/Li prescription of Bell 
et al. (2003) with a Chabrier (2003) IMF; 

Col. (19): logMbar, the baryonic mass of the galaxy in units of Mq . Section [s] discusses 
the sum of the stellar, Hi and gas masses to make up the total baryonic mass, 

Col. (20): logMx, the total mass of the galaxy computed at i? = -R23.5 according to 
Mt{R) = 2.33 X lO^RV^lJsin^ii) Mq where R is the radius in kpc and V^ot is the 
observed rotation velocity in kms~^. 

6. SDSS Comparisons 

Below, we compare our detailed structural parameters extracted from SDSS galaxy 
images with SDSS pipeline data products for the same galaxies. 

6.1. The SDSS Data Products 

The SDSS measurements of interest to us are the g, r, and i band Petrosian parameters. 
The Petrosian radius, i?p, is the radius at which ratio of the local surface brightness 
averaged over the annulus R^ is equal to 0.2 times the mean surface brightness within R^ as 
measured by the automated SDSS pipeline (Blanton et al. 2001; Strauss et al. 2002; Yasuda 
et al. 2001). The Petrosian magnitude, mp, is measured within the circular aperture of 
radius 2Rp. The Petrosian radii -Rpso and -Rpgo encompass 50% and 90% of the total light 
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measured within i?pQ 

We first find that the SDSS data products are burdened by the mis-identification 
of galaxy objects in the SDSS SQL query method, leading to erroneous measurements 
of radius and light for numerous objects. Some of the faulty measurements are due to 
the "shredding" caused by the deblending algorithm in the SDSS reduction pipeline. 
Overlapping objects are deblended to separate underlying components and extract proper 
measurements. Occasionally large galaxies may be interpreted as multiple systems and 
are thus "shredded" by the algorithm. The 1% "shredding" occurrence in the SDSS DRl 
reported by Blanton et al. (2001) is also found in DR7. We confirmed this "shredding" 
effect in the SDSS DR7; we found object-type mis-identifications for objects classified as 
galaxies that are either star-forming regions or foreground stars, both of which are more 
concentrated than a typical galaxy. 

To identify the deviant SDSS pipeline data, we plot in Figure [8] the logarithmic 
difference between the 90% and 50% light radii, i?p9o — -Rpso, in the r-band against both the 
magnitude mp in the upper panel (a) and the logarithmic Petrosian radius, Rp, in the lower 
panel (b). The clear separation at log-Rpgo — -Rpso = 0.55, shown as a dashed vertical line, 
is the clear distinction between correctly identified galaxies and the mis-identified compact 
objects. Objects located left of this line were discarded. 

We also examined the r-band Petrosian radii and their associated errors. Faulty 
computations of the Petrosian radius or magnitude are given a default -1000 error; those 
are shown as red crosses in Figure [8j Many failed calculation of Petrosian radius are given a 
SDSS default value of 3 arc seconds; hence the horizontal array of crosses near logi?p = 0.48 
in panel (b). These faulty measurements apply mostly to faint galaxies whose flux is below 



^The online SDSS archive tags are petroRad, petroR50, petroRpn and petroMag. 
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a nominal threshold within the Petrosian aperture. 

Thus, we weed out SDSS data based on the following criteria: i) mis-identified 
non-galaxian objects with log i?p9o — -Rpso < 0.55, and ii) failed computation of the Petrosian 
radius with errors of -Rp < arc seconds in the r-band. This eliminates 14% of the SDSS 
pipeline data products; we are left with a "clean" SDSS sample of 2605 galaxies. 

The computation of SDSS data products within a circular aperture compared to our 
structural parameters extracted from isophotal ellipse fitting is an additional source of 
confusion in the comparison of our data. Circular and isophotal apertures clearly measure 
different portions of the galaxy. Projection effects for a galaxy bulge and disk within a 
circular or elliptical aperture yield very different luminosity profiles (e.g., Bailin & Harris 
2009). 

Such inclination effects between the SDSS Petrosian half-light radius, -Rp5o, and 
the half-light radius, i?e; measured from isophotal fitting are explored in Figure |9} The 
Petrosian radius is systematically smaller than R^. More importantly, this offset increases 
with inclination, from red crosses {% < 50°) to blue circles {i > 75°) in Figure [9| as the 
Petrosian aperture progressively samples less light (compared to the isophotal aperture) 
at higher tilt. Slight curvature may also be detected in each inclination distribution but 
overall, there exists a direct mapping between -Rpso and -Re as a function of inclination. 
The scatter between a Petrosian radius and an isophotal radius is entirely dominated by 
inclination. 

6.2. Comparison of Radial Measurements 

We now compare the "clean" i-band SDSS Petrosian radii for 2605 galaxies with our 
suite of radial parameters, namely R^, -R23.5 and Rd computed for Sample A (Section [3]). 
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Figure 10 shows the raw Petrosian radial measures i?p, i?p5o and i?p9o from the SDSS 
pipehne against two of our own radial measurements -R23.5 and Rg, extracted from our 
2-band light profiles. The red dashed lines have slope unity. 

We first observe that inclination drives the scatter in all the correlations between 
Petrosian (i.e. circular) and isophotal (i.e. elliptical) radii although that effect is strongest 
for -RpSQ. Of course, no such dependence is seen for _Re versus -R23.5 as both measures are 
derived from isophotal fits. There is also no one-to-one correlation between our radial 
measures and the SDSS Petrosian radii. Therefore, any scaling relation based on SDSS or 
isophotal sizes ought to yield different parameters (slope, zero-point, scatter). 

The Petrosian radius Rp correlates well with the two other Petrosian radii -Rp5o and 
i?p90, as one might expect since the latter two derive from Rp. The bright end of the 
diagrams involving Petrosian radii could be slightly skewed to smaller values due to the 



reported over-sky subtraction errors for large objects (see Section 3.2). The effect of an 
over-subtracted sky on the measured isophotal radii, and fluxes, was demonstrated in 
Figure |3] for the sample galaxy UGC 5651. 

Table [9] presents the propagated error in radius AR^. and A_R23.5 as well as the 
magnitude differences Am23.5 and Amext due to sky uncertainty. A 0.2% sky error typically 
yields 1-2% radial scale variations or magnitude differences less than 0.03 mag. These 
small values are representative of our data; 1% sky errors would be dramatic but they are, 
fortunately, unrealistic. As observed earlier (Section 3.2), the isophotal parameters, -R23.5 



and m23.5 are least affected by typical sky background errors. 



Figure 10 highlights the comparison and differences between radii derived from circular 
and isophotal apertures. However, to decide which of these measurements yields the tightest 
galaxy relations requires an objective comparison against an uncorrelated variable, here 



chosen to be the deprojected rotational velocity, V. Figure 11 shows the distribution of 
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RV data for each i-band radial parameter; Sample A (all the available data) is shown with 
gray points overplotted whereas Sample D (best inclination and distance data) is in black. 
The top panels show RV scaling relations with isophotal radii and the disk scale length; 
similarly for the bottom panels with Petrosian radii for the clean SDSS sample. The red 
dashed line in each panel is an orthogonal linear fit, with bootstrap re-sampling (see C07), 
of the RV relations for Sample D. 

We evaluate the tightness of each fit with both the one-ajiv standard deviation about 
the best orthogonal fit and the Pearson r coefficient. The listed values of r and any are 
those for Sample D. 

The tightest RV relation is clearly that which involves -R23.5, in agreement with 
Saintonge et al. (2008) and SSll. We discuss their results in Section [7j Of the Petrosian 
radii, -Rpgo, yields the tightest RV relations. 



Figure 12 shows the variation of crjiy for the three gri radii shown in green, yellow and 
red colours respectively. The four sub-samples are plotted as A (circles), B (triangles), C 
(squares) and D (stars). The vertical dashed line separates our radial measurements from 
the SDSS Petrosian radii. 

Overall, Sample C and D yield comparably tight RV relations indicating that the 
inclination cuts in Sample D do not improve the trends for the RV relation significantly. 
Other RV relations with R^, R^ and the Petrosian radii are however improved when high 
and low inclinations are weeded out. Still, the biggest effect in obtaining a tighter RV 
relation results from eliminating galaxies with uncertain distances. 

Parameters derived from the extrapolated light profiles, such as -Re and R^, are 
especially sensitive to the vagaries of the "disk-like" exponential region in such profiles. 
Indeed both Re and R^ show a poorer scatter in Figure [l2j Rd is also affected by projection 
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effects, internal extinction by dust and stellar populations. For instance, there is a real 
dust extinction gradient from center (opaque) to edge (transparent) in galaxy disks. The 
measured scale lengths are larger, and central surface brightnesses are smaller than their 
intrinsic dust-free values. Empirical corrections (Giovanelli et al. 1995; Masters et al. 2003; 
Graham & Worley 2008) or corrections based on radiative transfer models (Popescu 
et al. 2005; Driver et al. 2007; Gadotti et al. 2010) for extinction effects on as a function 
of inclination have been proposed, but these remain grossly uncertain. 

The isophotal radius -R23.5 is the clear winner with this test, for all bands and 
all samples. The SDSS Petrosian radius, -Rpgo, also fares well though the unreliable 



object identification and deprojection effects for tilted galaxies (Section 6.1) make the 
interpretation of this parameter less simple than -R23.5- 

In seeking a direct match with galaxy formation models (e.g.. Mo, Mao, & White 1998; 
D07; Dll), numerous scaling relation studies have relied on as the fiducial galaxy size 
(e.g., C07 and references therein). This is because the scale length indicates the change 
of surface density with radius for any galaxy; as a global tracer of galaxy structure, Rd 
is a straightforward prediction of galaxy formation models. Different stellar populations 
however have different scale lengths. The space density of a given stellar population in a 
disk is mostly controlled by the angular momentum evolution of that disk (Foyle et al. 2008; 
Rov skar et al. 2008; Dutton 2009). R23.5, which measures the local surface density of 



a disk, is also sensitive to the local stellar populations (see Figure 13). The local space 



density is more sensitive to local star formation conditions than those averaged at -Rd- 

Whether -R23.5 or R^ is preferred for either the tightest RV relation that it yields or for 
a more direct connection with theory, the interpretation of the RV relation relies on a full 
understanding of the biases that each method entails. 
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6.3. Comparison of Luminosity Measurements 

In order to determine the best luminosity measure from circular (Petrosian) or elliptical 
apertures, we compare the clean SDSS Petrosian apparent magnitudes, mp, against our 



suite of derived apparent magnitudes m23.5, mgxt, nid and m2.2- Figure 14 shows these five 



2-band measurements against each other; the red dashed line has slope unity. 



The top row in Figure 14 shows that the best match with the Petrosian magnitude. 



mp, is obtained with m23.5. The top panels all show a departure from the one-to-one line 



at the bright end. This is still due to the "shredding" effect (see Section 6.1) for large 
galaxies. For non-shredded systems, the Petrosian fiux accounts for ~ 98% of the total fiux 
of a galaxy (see Shen et al. 2003). 

The disk magnitude mj shows the largest scatter with other luminosity measurements, 
as it is tied to the uncertain disk scale length via Equation ^ and does not account 
for the bulge luminosity. A superb correlation between m23.5 and m-ext exists; m2.2 however 
yields noisier correlations. 



In Figure 15 we examine which i-band luminosity parameter yields the tightest LV 
scaling relation. The dashed red line is the orthogonal bootstrap fit through Sample D, 
shown with black points. The full Sample A is shown in gray. The extrapolated luminosity, 
Lext; and the isophotal luminosity, -^23.5, yield mathematically tightest LV relations; 
however, the Petrosian radius, Lp, and L2.2 are practically as good. 



Figure 16 shows the la standard deviation for the LV fits in each gri band and for 
each galaxy sub-sample A-D. Sample D at i-band yields tightest relations, as does the use 

of Lext- 

Overall, all the luminosity metrics, with the exception of the disk luminosity, Lj, 
yield comparably tight LV relations. Technically the total extrapolated luminosity of the 
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galaxy Lcxt and L23.5 yield slightly tighter LV relations than our other tested luminosities. 
Luminosity measurements are rather stable due to their global (cumulative), rather than 
local, nature. 

We have seen in Table [o] the effect of a ±0.2% and ±1% sky uncertainty on the radial 
and apparent magnitude measurements for the galaxy UGC 5651. The luminosity within 
2.2 disk scale lengths, m2.2, is mostly affected by uncertain disk scale lengths. Extrapolated 
magnitudes are less sensitive with ~ 0.03 mag uncertainty for a typical 0.2% sky error. 
Despite small luminosity errors, we stress the importance of well-measured sky levels, and 
caution against SDSS luminosity estimates for the brightest and biggest galaxies in our 
sample. 



7. VRL Relation 

We now use the best galaxy structural parameters, -R23.5, -Rd, -^cxt and L23.5 as 
determined in Section |6| and the corrected galaxy rotational velocity VJ-ot (hereafter V) for 
a detailed analysis of VRL scaling relations. 

We construct various scaling relations VL, RL and VR which we refer to collectively as 
the VRL relation. We also adopt orthogonal linear fits to model our scaling relations. Fit 
differences result from different modeling techniques; e.g., bisector fits yield a steeper slope 
than orthogonal fits. Given that the VRL scaling parameters are not fully independent 
(e.g., correlated via inclination and distance), neither a forward or inverse fit would work. 
Bisector fits, which average the forward and inverse fits, are therefore also incorrect by 
construction. Indeed, the bisector fit of a perfectly uncorrelated distribution of two variables 
has an absolute slope of one (C07; Hogg et al. 2010). Lacking an accurate covariance 
matrix, we adopt the orthogonal fit, as did C07 and SSll, as the least biased of the suite 
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of numerical fitting methods. The VRL parameter errors discussed in Section 4.1 are used 
for, but play little role in, our orthogonal fits. 



Figure 18 shows the logarithmic form of the VRL relation, using the scaling parameters 
-^ext; -R23.5 and V, for our moderate inclination, best determined distance Sample D in 
black. The full Sample A is shown in gray. The red line is the orthogonal fit to Sample D's 
parameter distribution and the 2-a bootstrap errors are shown as dotted lines. See C07 for 
details about the fitting procedure. The direction and magnitude of a 20% distance error 



are shown with an arrow in each panel's corner (see Section 7.4 ) 



Tables [9j - |9j present the fit slopes and zero points, as well as the forward Y\X and 
inverse X\Y 1-a deviations and overall Pearson r coefficients of the gri-ha.nd VRL relations. 
The crx\Y and <7y\x standard deviations are presented to facilitate comparisons with other 
authors and/or theoretical models (D07; Dll). Any observational limit (e.g., magnitude 
cut) will bias the estimate of (Jx\y against aY\x- The X\Y and Y\X sigmas also differ in 
proportion to the slope of each scaling relation. 

To compare our standard deviations with those of other authors, let us write the 
slope and scatter transformations for systems relying on magnitude or luminosity. Ignoring 
zero-points, the VL relation (logV = fey^logL) is directly related to the VMag relation 
{logV = bvMagMag) via Mag = — 2.51ogL or byi = —2.5 x by Mag- Similarly, the scatter 

<^VL = <^MagV I^MagV- 

Given that previous scaling relation studies (e.g., C07; AR08; D07; SSll) have used 
both i?d and -R23.5 as radial metrics, the various scaling rleations that we present in Tables 
|9]-[9]will showcase those two radii where relevant. 
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7.1. The V-L (Tully-Fisher) Relation 



The VL relations from Section pl3\ are modeled as logl^ = a + blogL. We tabulate in 
Table [9] our best fit values for a and b in the i-band as a function of luminosity (whether 
Lcxt or -L23.5). For our best i-band Sample D with 652 galaxies, we find 

yocL°-™ V^ocL°3™. (15) 



The slopes and scatters for the multi-band gri VL relations are tabulated in Table [9j 
The tightest relation is indeed obtained in the i-band, and slopes are slightly shallower from 
g to i band. 

These i-band slopes are identical to those reported by M06, C07 and SSll (among 
others). The consistency of the VL slopes from these and other studies based on rather 
different samples attests to the robustness of the VL relation and its independence to 
numerous "third parameters" (e.g., Courteau & Rix 1999; D07). The VL scatters reported 
in the literature may differ more significantly as a result of various selection and/or 
"modeler" biases. Since we compare our scatters with those of M06, COT and SSll below, 
a brief reminder about each sample is warranted here: 

C07 studied the VRL scaling relation of 1300 late-type galaxies with a mixture L23.5 
and Lext luminosities and Ha and Hi rotation curve velocities. M06 used a subset of 807 
so-called "template" cluster galaxies from the S07 catalog. They extracted their own /-band 
total magnitudes truncated at 8-Rd and corrected for Galactic and internal extinction 
and k-cosmological term. These magnitudes were also corrected for morphologjj^ as well 
as incompleteness bias. SSll used the entire S07 catalog (template and non-template 



^While appropriate for distance scale studies, a morphological correction to galaxy mag- 
nitudes removes all astrophysical signatures of morphology on scaling relations. 



- 30 - 



galaxies). Like C07, they did not correct for morphological type and incompleteness bias, 
but unlike C07, their final fits invoked outlier clipping. The Hi line widths for M06, SSll 
and the current study all come from the same source (S07). 

In order to compare scaling relation scatters obtained from orthogonal fits, we were 
fortunate to use the original data files from M06 (publicly available; 807 galaxies) and SSll 
(A. Saintonge, private comm.; restricted set of 665 template galaxies). 

For our sample, we find cTy^" = 0.074 (Table |9]). Meanwhile, our orthogonal VL fits 
to the M06 and SSll linewidths and /-band magnitudes of their template galaxies yield 
ay^° = 0.065 in both cases. Considering the similar samples, we attribute most of the 
scatter difference with M06 to their morphological correction and with SSll to their outlier 
clipping procedure. Our tests have also revealed scatter differences of 0.02 dex between our 
respective orthogonal fitting engines (i.e., same data, different software). 

CO 7 found an even smaller observed scatter of avL = 0.057. Note that CO 7 did not 
correct for morphological differences or for incompleteness bias, or sigma-clip their scaling 
relations. We speculate that the scatter difference is here due to the morphological make-up 
between the S05/S07 samples and the ones collected in C07. The former is morphologically 
broader than C07. For instance, restricting our sample to only Sc galaxies reduces the VL 
scatter by 0.005 dex. Tailoring one's data set according to specific criteria will obviously 
affect the final scaling relations and complicate comparisons amongst even similar studies. 

Other related results are also found in the literature. Kannappan et al. (2002) obtained 
a 5-band VL slope of 0.29 and scatter of ayi = 0.08 for 68 spiral galaxies. AR08 found 
a i?-band VL slope of 0.314 ± 0.015 and an intrinsic scatter of avL = 0.063 for 76 disk 
galaxies. The TF study of 162 disk galaxies with Ha rotation curves and SDSS imaging 
by Pizagno et al. (2007; hereafter P07) allows a more direct comparison with our SDSS 
based results. Their derived bivariate fit gri VL slopes are bvL,g = 0.340, bvL,j- = 0.338 and 
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bvL,i = 0.325, all steeper than ours (Table |9]). The scatter in each relation (JvL,g = 0.073, 
o'vL,r = 0.063 and crvL,i = 0.061, is also smaller than ours, partly due to their steeper slopes, 
different fitting methods (bisector vs orthogonal), and their more careful sample selection. 

While many of the studies above either seem to agree or differ slightly, we stress 
that the comparison of scaling relation parameters depends highly on the selection and 
manipulation of each data sample. 



7.2. The R-L (Size-Luminosity) Relation 

Much like the VL relation, the size-luminosity RL relation is a fundamental constraint 
to galaxy formation models (P07, AR08, D07, Dll). We compute the RL relation for the 
combinations of our radial parameters Re, R23.5 and R^ with the luminosity parameters 



Lcxt) -^23.5; and L2.2- The 1-cr scatter of each RL relation is shown in Figure 17, where 
each subsample A-D is separated by point type and each band by colour (see caption). Of 
our three radial measurements, the isophotal radius -R23.5 yields the tightest RL relation. 

We present the four combinations of the RL relation with R23.5, Rd, -^23.5 and L^^t in 
Table |9} We have also computed fits using the corrected radii -R23.5 R^, as decribed 
in Section |4j however, these corrections provide no gain in scatter reduction. In view of 
their uncertain nature, we use uncorrected radii in what follows. We express our results as 
R (X L^ for our best i-band Sample D of galaxies with uncorrected radii, 

D „ r0.41±0.02 p „ r0.36±0.20 (^p_\ 

oc Lp^t Rd (X L23.5 (16) 

p „ r0.44±0.07 p ^ r0.41±0.07 /i 



As both radii and luminosities depend on the distance estimate (Section |4j), distance errors 
contribute weakly to the scatter in the RL relation over all subsamples A-D. We show in 



Figure HiSl and in Section 7.4 that distance errors scatter along the direction of the RL 
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relation, thus introducing minimal spread to the relation. The larger Sample B (moderate 
inclination galaxies only) is thus often just as tight as the smaller, distance-pruned. Sample 
D of galaxies. 



The measured slopes in Equation (17) are in broad agreement with previously 
reported values (e.g., de Jong & Lacey 2000; Shen et al. 2003; C07; SSll). The gri band 
dependence of the RL relation is shown in Table [9] As in the VL relation, the RL relation 
slope decreases from bluer to redder bands. Bandpass effects are least dominant for the 
combination of i?d and ^^23 5 and largest for the pair -R235 ^"^^ -^Lt- 

Use of scale lengths to construct the RL relation yields weak correlations with a 
Pearson coefficient r ~ 0.7 compared to r ~ 0.9 for isophotal radii. The larger scatter 
of the scaling relations based on R^ partially stems from the larger uncertainty in the 



measurement of the disk scale length itself (Section 6.2). SSll found a similar result; their 
isophotal radii measurements yield an RL r ~ 0.86 while C07's RL relation based on 
disk scale lengths has r ~ 0.65. Our use of -R23.5 is thus an improvement over C07 and a 
confirmation of SSll's results (see also C96). 

The different scatters for our RL relations are listed in Table [9| with the tightest value, 
o'rl = 0.068, obtained for the combination of -R23.5 and L^^^. The RL scatter ajiL = 0.05 
reported by SSll for their isophotal RL relation is tighter than our or any other previous 
similar assessment (e.g., Shen et al. 2003; Pizagno et al. 2005; C07; AR08; see Fig. 3 of 
SSll). Eliminating SSll's outlier clipping procedure increases their reported RL scatter to 
crl = 0.06, in closer agreement with our value. 
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7.3. The R-V (Size- Velocity) Relation 

The third scaling relation combines the galaxy size, R, and velocity, V. Table [9] shows 
our orthogonal fit parameters for the RV relation. For instance, the i-band RV slopes for 
Sample D are: 

i?d OC V^l-82±0.14 ^^^^ ^ ^1.52±0.07_ ^^g) 

The observed log scatter of the -R23.5 — V relation is 0.154 dex. Table |9] gives the slopes 
and scatter for Sample D in all three gri bands. 

We can compare the i-band RV parameters with those of SOS and SSll. SOS's 
bivariate fit to a sub-sample of 699 S07 galaxies yielded i?23 5''^'^'^'^ oc V^i-24±o.oi ^^^^ ^ scatter 
'^w'^''^'^'^ = 0.07. In order to compare directly with SOS, we compute bivariate fits for the 
RV relation of our Sample D yielding: Rfi""^ oc •l/i-37±o.oi ^bisec ^ 0^45 Qur and SOS's 
results are significantly different. Comparison with SSll's similar, yet larger, dataset yields 
somewhat different results: i?23.5^'°'^*'^° Vi-44±o.o2 ^ scatter a^Y°'^^^° = 0.123. SSll's 
RV slope is statistically equivalent to ours; the smaller scatter results again from outlier 
clipping. 

Due to the larger scatter in their RV relation, C07 derived their RV slope by combining 
the VR and RL relations and requiring that if oc and R oc then R oc The 
C07 RV relation had a (log) slope of 1.10 ± 0.12 dex. Our RV relations listed in Table [9] 
are all steeper and obey intrinsically the rule R oc V"'^^ (with R oc V^-^^Z-^''), this sense, 
this study supercedes C07. 

i?23.5 remains the ideal radial metric for the tightest correlation with V but, of the 
three scaling relation combinations {VL, RL and RV), the RV relation is broadest. This 
weakens its appeal for cosmological (e.g., distance-measuring) applications. 
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7.4. Distance Dependence of the VRL Relation 

We calculate the effect of a percent distance error AD/D in the LV, RL and RV 
planes as, 

AlogL=^^ = 0.868^ (19) 

^ 2.51n(10)D D ^ ' 

A\ogR = , = 0.434—. (20) 

^ ln(10)D D ^ ' 



We show the 20% distance uncertainty in the bottom right corner of each Figure 18 
plot. Note that the propagated radial uncertainty A log R is twice as large as the uncertainty 
in luminosity A log L and that distance uncertainties in the RL relation move points along 
a slope of 1/2. Since the slope of the RL relation is ~ 0.5 (Table |9]), distance uncertainties 
play no role in the RL relation scatter. 

Recall that the distances used to calculate luminosity and physical radius of our galaxies 



are derived from the CMB redshift velocities listed by SOS and S07 (see Equation (11)) 



SOS and S07 also present "peculiar velocity corrected" distances (in Equation (11), 
K) = VcMB + Vpec); these were computed to minimize the scatter in the TuUy-Fisher (VL) 
relation for each galaxy cluster in their sample. The peculiar velocity estimates can thus be 
very large and even unrealistic as they absorb all possible sources of scatter. We therefore 
caution against the use of such velocities in scaling relation studies. 



Figure 19 shows the VRL relation computed with both Vcmb distance data (black 



points and lines) and peculiar velocity data (red points and lines). Standard statistics 



are shown. The Vcmb data in black are the same as in Figure 18, Note the significant 
difference in slopes due to an over-interpretation of the scatter as being due solely by 
peculiar velocities. 



- 35 - 



7.5. Dependence on Morphology, Colour and Surface Brightness 

Much effort has been invested in the study of the VRL relation scatter, with particular 
emphasis on the Tully-Fisher VL relation. The search for the "third" parameter has 
included, among others, morphology, colour, mass-to-light ratios and surface brightness 
(Zwaan et al. 1995; CR99; Verheijen 2001; Kannappan et al. 2002; C07; SSll to cite a few). 
In this section we highlight the VRL dependence as a function of morphology, colour and 
surface brightness. The full quantitative dependence analysis will be presented elsewhere. 



Figure |20] presents the morphological dependence of the VRL relation for Sample A; 
the RL and RV relations on the left and right figures use -R23.5 and R^, respectively. The RL 
slope are expected to become shallower from early to late-type galaxies (Shen et al. 2003), 
as we see for R^. The isophotal RL relation is however free from such a trend. Fitting 
different VRL relations for different morphologies is beyond the scope of this work; however 
we do confirm the same trends as measured by C07. 

Galaxy colour has been examined as a possible agent in the VRL scatter, with larger 
variations in bluer bands due to star formation (see Kannappan et al. 2002; C07; POT). We 



show the colour dependence of the VRL relation in Figure 21 for Sample D only for clarity. 



We observe a slightly larger scatter of the VRL relation for bluer galaxies. Clearly the full 



Sample A (shown in Figure 20 ) shows the largest scatter in the redder, early-type galaxies. 



Overall this confirms that redder, brighter galaxies are faster rotators and that the 
largest galaxies are the brightest and reddest late-type galaxies. 



Figure 22 shows the dependence of the VRL relations on surface brightness. The 
effective surface brightness, /ig, is colour-coded with red for highest surface brightness 
(HSB) galaxies and blue for low surface brightness (LSB) galaxies. The RL relation 
illuminates the conjunction of high concentration galaxies (small R, large L) with HSBs 
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and low concentrations with LSBs. One might naively expect that HSBs are faster rotators 
than more extended LSB galaxies for a given luminosity, though no dependence of the VL 
relation on surface brightness has ever been found (e.g., Sprayberry et al. 1995; Zwaan 
et al. 1995; CR99; D07; SSll). Note that if the velocity were measured at relatively small 
radii, a surface brightness dependence of the VL relation would be expected (e.g., Catinella 
et al. 2007). VL and RL residuals are only weakly correlated indicating that for a given 
luminosity, the disk size does not affect the rotation speed. The intrinsic scatter in the 
VL relation may then be attributed to scatter in the dark matter fraction and stellar 
mass-to-hght ratio (D07). 

The non-correlation of in the VL residuals has often been interpreted as proof for a 
correlation of surface brightness with dark matter content. The latter is however difficult to 
assess given large uncertainties in stellar population and galaxy formation models. Various 
factors contributing to the VRL scatter include star formation, the initial mass function of 
stars, halo and disk spin parameters, feedback, adiabatic contraction of the halo, correlation 
of halo parameters and more (Navarro & Steinmetz 2000; Firmani & Avila-Reese 2000; D07; 
AR08; Sll, Dll). A detailed, quantitative model analysis of these data will be presented 
elsewhere. 



8. The Bciryonic TuUy-Fisher Relation 

The Baryonic TuUy-Fisher relation (BTF) relates the total baryonic (stellar + gaseous) 
mass of a spiral galaxy to its rotational velocity. For their intermediate-size samples 
(a few hundred galaxies), McGaugh et al. (2000), Trachternach et al. (2008), and Stark 
et al. (2009) found that the BTF is a more "fundamental" relation than the luminous 
(TuUy-Fisher) or stellar mass TuUy-Fisher (STF) relation. 
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We can test these postulates with SOS's larger compilation since it includes gas masses 
for thousands of galaxies. We can indeed sum the stellar masses, computed from SDSS 
luminosities and stellar population models (e.g., Bell & de Jong 2001), with the SOS gas 
masses to obtain a baryonic mass for each galaxy. 

McGaugh et al. (2000) also noted a kink in the luminous TF relation of the low-mass 
galaxies with V ^ 100 kms~^; those would have systematically lower luminosities than 
predicted by the TF fit at higher mass. By considering the total baryonic mass instead, 
they showed that this effect disappears and galaxies of both low and high baryonic assume 
the same BTF relation. This result is not confirmed by Geha et al. (2006) who do not see 
a deviation from the luminous high mass TF relation for a sample of 101 very low-mass 
dwarf galaxies all with V ^ 100 kms~^. Our sample does not permit a critical assessment 
of the McGaugh versus Geha dilemma as it includes too few galaxies with V ^ 100 kms~^ 
and Mgas > MF] However, it is worth computing the scatter and slope of the BTF down to 



^gas _ 

those velocities. To our knowledge, ours is the largest BTF sample to date. 

For the computation of our galaxies' STF and BTF relations, we use as the fiducial 
luminosity parameter our extrapolated light measurement L^^^ (Section |6.3 ) which provides 



the tightest correlation with rotational velocity (Section [T]). We convert our {g — i)ext galaxy 
colours into an i-band stellar mass-to-light ratios, M^,/Li, using the prescription of Bell 
et al. (2003) with a -0.10 offset to reproduce the Chabrier (2003) IMF. 

The neutral gas mass is calculated from the corrected Hi line flux measurement, Sc, in 
SOS following Haynes & Giovanelli (1984): 

Mhi = 2.3SS X lO^ScD^ (Mq), (21) 



^Note that the tension with Geha et al. (2006) is reduced if the improved data base of 
Stark et al. (2009) is used instead of McGaugh et al. (2000). 
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where D is the distance of the galaxy in kpc (see Equation (11)). The contribution from 



hehum and other metals to the gas mass is approximated as in McGaugh et al. (2000); 

Mgas = 1.4Mhi (Mq). (22) 
The baryonic mass is the sum of the stellar and gas contributions, 

Mbar = + Mgas (Mq) (23) 

The baryonic sample contains ~5% fewer galaxies than the STF and TF sample, given 



the limited availability of Hi flux measurements. Figure 23 shows the STF (top left) and 



BTF (top right) relations, respectively; the orthogonal fit results are tabulated in Table |9 



for our galaxy Sample D. We show in the left and right bottom panels of Figure 23 and in 
the last two entries of Table [9] the gas contributions to the total stellar and baryonic masses 
expressed as log(Mg/M*), the logarithmic ratio of gas to stellar mass. For our best Sample 
D, the STF and BTF relations are fit as: 

V oc mO-25±o-oi V cx M^f (24) 

with measured scatters ayu* = 0.072 and o-vM^ar ~ 0.076, respectively. 



These scatters are comparable to those of the basic VL relation (Section 7.1), which is 
remarkable given the greater number structural parameters involved. Given the uncertainty 
in assessing the assumed errors on M^,/L, L, one would be hard-pressed to identify any 
differences between the STF and BTF scatters; they are both comparable. However, while 
the observed scatters of the TF, STF and BTF relations are comparable, both the STF and 
the BTF relations may have intrinsically smaller scatter than the VL relation and might 
thus be more fundamental. 

Figure 24 shows the dependence of the STF and BTF relations of the {g — z)ext colour, 
as well as their gas contributions, for sample D. The bottom panels show that the ratio of 
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gas to stellar and baryonic mass is well-defined by colour: larger more massive (luminous) 
galaxies have a smaller gas fraction than smaller, bluer galaxies. The gradient in gas mass 
fraction, from the least to most massive galaxies, causes most of the steepening of the BTF 
relation compared to the STF. 

The Pearson r correlation coefficient suggests a mildly tighter STF relation (versus 
the BTF). This may be partially attributed to the distance-dependence of the Hi gas 



mass calculated in Equation (21) which introduces additional scatter into the BTF. It is 
also interesting to note that for their semi-analytic model. Button, van den Bosch, & Dekel 
(2010) showed that galaxy star formation rates (SFR) are more tightly correlated with 
stellar mass than baryonic mass (see their Fig. 15). However, to determine which of the STF 
or BTF relation is more fundamental will require far more accurate measurement errors for 
the stellar and gas masses. For instance, our gas mass estimates include a correction for 
helium and heavier elements but do not include molecular hydrogen or the warm ionized 
gas which could increase the baryonic fraction and steepen the BTF slope slightly (see 
Fukugita & Peebles 2004; Bregman 2009; McGaugh et al. 2010). A discussion of these 
effects is beyond the scope of the present study, but should be addressed elsewhere soon. 

The uncertainty in the stellar mass-to-light ratio, M^,/L, propagates through both STF 
and BTF relations and different calibrations will therefore yield different slopes and scatters. 
McGaugh et al. (2005) explored three variants of M^/L and their effect on the scatter of 
the resulting BTF. Their first M^^/L assumed a maximum disk hypothesis. The second 
method used the B — V colour prescription of Bell et al. (2003), and the third estimate 
was based on modified Newtonian dynamics (MOND). McGaugh et al. (2005) concluded 
that the MOND M*/L's yield the tightest (r = 0.99) BTF with a slope of 0.25. However, a 
^bar ~ Vfiat relation with slope 0.25 and zero scatter is intrinsic to MOND, and its recovery 
via MOND M*/L's should not be surprising. Furthemore, the MOND prediction applies for 
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McGaugh's velocity measure Vfiat (McGaugh 2000), whereas our BTF relation relies on Hi 
line widths, W^^q. In dwarf galaxies, Wf^q is likely to under-estimate Vfiat, while in massive 
spirals FFpso is likely to over-estimate Vfiat (Verheijen 2001). The BTF slope based on Hi 
line widths should thus be steeper than that based on Vfiat measurements, as observed. 
McGaugh (2005) certainly found a range of BTF slopes from 0.25 to 0.33 depending on the 
velocity estimator. 

Our computation of the BTF relation used the sub-maximal M^/L values of Bell 
et al. (2003). Adopting a maximal disk decreases the BTF slope from 0.292 to 0.290; the 
effect is thus neghgible. Bell & de Jong (2001) find the same slope of ~0.29 for the BTF 
derived with maximal M*/L's in the B, R, I and K bands. 

AR08 also studied 76 disk galaxies with masses derived from re-scaled M^/Us from 
Bell & de Jong (2001) and obtained a STF slope of 0.274±0.012 and a BTF slope of 
0.306±0.012, in good agreement with our own results. A shallower slope is however found 
by Gurovich et al. (2010) who fit a bivariate BTF for a small sample of 21 galaxies with 
Hi line widths and the V- and i?-band luminosities to find a BTF slope ~ 0.26. Here is a 
case where the use of fine widths does not yield a shallow BTF slope; the small size and 
bivariate fitting method however complicate the interpretation of this different result. To 
our knowledge, use of the velocity measure Vfiat always yields shallower BTF relations. 

A more extensive census of BTF relation slopes is presented in Gurovich et al. Overall, 
BTF studies report V — M^ar slopes in the range 0.25-0.33. However, care in comparing 
BTF relations from different authors, using different definitions of Vdrc (e.g., whether Ymax, 
W-F50, Vfiat, etc.) and different fitting methods (e.g., bivariate vs orthogonal fits) must 
always be taken. 
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9. Conclusion 



We have compiled a catalog of 3041 spiral galaxies within cz < 28000 kms~^ with 
multiband photometric parameters extracted from SDSS images and corrected Hi velocity 
widths and distances taken from Springob et al. (2005, 2007). We have extracted g^r and 
i-band surface brightness profiles from SDSS images and derived well-cahbrated radial and 
luminous measurements from each galaxy. Our main findings are as follows: 

1. Our radial and luminosity parameters extracted from isophotal fits of SDSS images 
yield a scatter improvement of the VL and RL relations of ~ 8% and ~ 30% compared 
to similar relations constructed with SDSS DR7 Petrosian parameters. 

2. We find that the i-band isophotal radius -R23.5 correlates best with the galaxian 
rotational velocity, V^ot- While i?23.5 is a more robust measurement than the scale 
length i?d, for instance in terms of the RL relation scatter, R^ holds greater 
significance for comparisons with galaxy formation models. 

3. Both the extrapolated and isophotal i-band luminosities provide the tightest 
correlations with galaxian rotational velocity. The luminosity within 2.2 stellar disk 
scale lengths (peak of the baryonic rotation curve) is a poorer tracer of T^ot- 

4. We construct the VRL relation in the gri bands finding general agreement with the 
slopes and scatters reported elsewhere. Our main scaling relations in the i-band are 
summarized as: 



V (xL\ 



0.27±0.01 
'23.5 



(25) 



-R23.5 L. 



0.40±0.01 
'23.5 



i?d (X L. 



0.35±0.02 
'23.5 



(26) 



-R23.5 oc V 



1.52±0.07 



i?d ocl/ 



1.82±0.14 



(27) 



(28) 
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We find little correlation of the VRL scatters with morphology, colour or surface 
brightness. Scahng relations for the SDSS g and r bands are provided in the text. 

5. We have transformed our galaxy colours into stellar masses and added these to the gas 
mass from the available Hi fluxes. These masses yield the largest stellar TuUy-Fisher 
(STF) and baryonic TF relations (BTF) to date. For the STF and BTF, we find 

V (X M°-2^±°-°i V oc (29) 

both with high Pearson r correlation coefficients of ~0.9. 

The observed scatters of the TF, STF and BTF relations are all comparable. To 
decide which of the STF or the BTF relation is more fundamental requires a detailed 
error investigation which is beyond the scope of the present paper. Our reported BTF 
slope (0.29) matches that of previous similar investigations based on Hi or Ha hne 
widths. Use of the velocity measure Vfiat would yield a shallower BTF relation slope 
of 0.25. 
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Table 2. Photometric Uncertainties for UGC 5651 with ±0.2% and ±1.0% Sky Error 



Radius (%) 


+0.2% 


-0.2% 


+1.0% 


-1.0% 


m (mag) 


+0.2% 


-0.2% 


+1.0% 


-1.0% 




-1.70 


1.88 


-5.62 


21.00 


Amext 


0.03 


-0.03 


0.15 


-0.16 




-0.69 


1.00 


-3.93 


5.14 


A m23.5 


0.01 


-0.01 


0.07 


-0.08 




-9.65 


6.19 


-24.26 


27.16 


Amd 


-0.03 


0.00 


-0.04 


-0.09 












A m2.2 


0.08 


-0.05 


0.21 


-0.20 
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Table 3. Orthogonal Fits for the i band VL Relation 
log V — a + b log L 



Sample 




a± Aa 


b± Ab 






Pearson r 




















A 


3041 


-0.688 ± 


0.060 


0.279 ± 


0.006 


0.141 


0.504 


0.711 


B 


1725 


-0.829 ± 


0.079 


0.292 ± 


0.008 


0.136 


0.467 


0.729 


C 


1076 


-0.670 ± 


0.054 


0.279 ± 


0.005 


0.075 


0.269 


0.885 


D 


652 


-0.761 ± 


0.086 


0.287 ± 


0.008 


0.075 


0.261 


0.872 


T/ 1 

►^rot 1-^23.5 


















A 


3041 


-0.503 ± 


0.056 


0.263 ± 


0.005 


0.139 


0.529 


0.718 


B 


1725 


-0.619 ± 


0.074 


0.274 ± 


0.007 


0.134 


0.491 


0.737 


C 


1076 


-0.525 ± 


0.050 


0.266 ± 


0.005 


0.074 


0.278 


0.888 


D 


652 


-0.595 ± 


0.079 


0.273 ± 


0.008 


0.074 


0.269 


0.877 
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Table 4. Orthogonal Fits for the gri band VL Relations of Sample D 

log V — a + b log L 



Band 


N 


a± Aa 


b± Ab 






Pearson r 
















g 


652 


-0.995 ± 0.108 


0.312 ± 0.010 


0.083 


0.266 


0.844 


r 


652 


-0.808 ± 0.092 


0.294 ± 0.009 


0.077 


0.264 


0.864 


i 


652 


-0.761 ± 0.086 


0.287 ± 0.008 


0.075 


0.261 


0.872 


T/ 1 

rot 1-^^23.5 














g 


652 


-0.691 ± 0.095 


0.286 ± 0.009 


0.082 


0.287 


0.848 


r 


652 


-0.624 ± 0.083 


0.278 ± 0.008 


0.075 


0.271 


0.871 


i 


652 


-0.595 ± 0.079 


0.273 ± 0.008 


0.074 


0.269 


0.877 
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Table 5. Orthogonal Fits for the i band RL Relation with Uncorrected Radii 

log R = a + b log L 



Sample N a± Aa fe± Ab grl (Tlr Pearson r 





A 


3041 


-3 


312 


± 





075 


0.377 


± 





007 





171 





454 





743 




B 


1725 


-3 


531 


± 





104 


0.396 


± 





010 





165 





417 





763 




C 


1076 


-3 


028 


± 





118 


0.349 


± 





oil 





153 





440 





744 




D 


652 


-3 


454 


± 





179 


0.387 


± 





017 





151 





392 





751 


Rd 


-^23.5 
A 


3041 


-2 


821 


± 





083 


0.332 


± 





008 





185 





557 





689 




B 


1725 


-2 


974 


± 





120 


0.345 


± 





012 





182 





527 





702 




C 


1076 


-2 


677 


± 





121 


0.317 


± 





012 





162 





513 





704 




D 


652 


-2 


996 


± 





193 


0.345 


± 





019 





163 





475 





700 


R23.5 


-^ext 

A 


3041 


-3 


251 


± 





041 


0.419 


± 





004 





096 





230 





918 




B 


1725 


-3 


500 


± 





041 


0.441 


± 





004 





073 





166 





953 




C 


1076 


-2 


925 


± 





060 


0.390 


± 





006 





089 





228 





914 




D 


652 


-3 


345 


± 





067 


0.427 


± 





006 





071 





165 





944 


R23.5 


-^23.5 

A 


3041 


-2 


948 


± 





039 


0.393 


± 





004 





098 





249 





915 




B 


1725 


-3 


147 


± 





044 


0.410 


± 





004 





079 





191 





946 




C 


1076 


-2 


692 


± 





058 


0.370 


± 





006 





091 





247 





908 




D 


652 


-3 


052 


± 





075 


0.402 


± 





007 





076 





189 





935 
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Table 6. Orthogonal Fits for the gri band RL Relations of Sample D 

log R — a + b log L 



Band N a± Aa 6± Ab grl (Tlr Pearson r 



Ra 


l-^ext 
g 


652 


-3.831 


± 





183 


0.424 


± 


0.018 





145 





343 


0.763 




r 


652 


-3.476 


± 





169 


0.389 


± 


0.016 





144 





371 


0.761 




i 


652 


-3.454 


± 





179 


0.387 


± 


0.017 





151 





392 


0.751 




-^23.5 
g 


652 


-3.037 


± 





191 


0.351 


± 


0.019 





165 





470 


0.674 




r 


652 


-2.963 


± 





177 


0.342 


± 


0.017 





158 





464 


0.701 




i 


652 


-2.996 


± 





193 


0.345 


± 


0.019 





163 





475 


0.700 


-R23.5 


1 

l-^ext 
g 


652 


-3.917 


± 





072 


0.474 


± 


0.007 





071 





150 


0.945 




r 


652 


-3.471 


± 





066 


0.438 


± 


0.006 





069 





157 


0.947 




i 


652 


-3.345 


± 





067 


0.427 


± 


0.006 





071 





165 


0.944 


-R23.5 


-^23.5 
g 


652 


-3.446 


± 





078 


0.434 


± 


0.008 





073 





167 


0.942 




r 


652 


-3.143 


± 





073 


0.410 


± 


0.007 





074 





181 


0.937 




i 


652 


-3.052 


± 





075 


0.402 


± 


0.007 





076 





189 


0.935 
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Table 7. Orthogonal Fits for the i band RV Relation 
logR^a + blogV 



Sample 


N 


a± Aa 


b± Ab 






Pearson r 




















A 


3041 


-2.944 ± 


0.161 


1.606 ± 


0.073 


0.312 


0.194 


0.428 


B 


1725 


-2.972 ± 


0.212 


1.614 ± 


0.096 


0.310 


0.192 


0.437 


C 


1076 


-3.087 ± 


0.196 


1.661 ± 


0.089 


0.226 


0.136 


0.592 


D 


652 


-3.484 ± 


0.315 


1.821 ± 


0.141 


0.240 


0.132 


0.563 


R23.5 1 Kot 


















A 


3041 


-1.901 ± 


0.094 


1.356 ± 


0.042 


0.227 


0.168 


0.610 


B 


1725 


-1.871 ± 


0.109 


1.339 ± 


0.049 


0.212 


0.158 


0.655 


C 


1076 


-2.151 ± 


0.106 


1.475 ± 


0.048 


0.158 


0.107 


0.761 


D 


652 


-2.287 ± 


0.149 


1.519 ± 


0.066 


0.154 


0.101 


0.762 
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Table 8. Orthogonal Fits for the gri band RV Relations of Sample D 

logR^a + blogV 



Band 


N 


a± Aa 


b± Ab 


ctrv 




Pearson r 


-n-d 1 I' rot 














g 


652 


-3.303 ± 0.310 


1.733 ± 0.139 


0.229 


0.133 


0.565 


r 


652 


-3.303 ± 0.293 


1.732 ± 0.132 


0.230 


0.133 


0.563 


i 


652 


-3.484 ± 0.315 


1.821 ± 0.141 


0.240 


0.132 


0.563 
















g 


652 


-2.439 ± 0.155 


1.539 ± 0.069 


0.162 


0.105 


0.742 


r 


652 


-2.313 ± 0.148 


1.514 ± 0.066 


0.155 


0.102 


0.758 


i 


652 


-2.287 ± 0.149 


1.519 ± 0.066 


0.154 


0.101 


0.762 



Table 9. Orthogonal Fits for the Stellar Mass and Baryonic TuUy-Fisher Relations 

Y\X, logr ^a + b\ogX 



Sample D 


N 


a ± Aa 


6± A6 






Pearson r 


V\M, 


652 


-0.429 ± 0.068 


0.253 ± 0.006 


0.072 


0.292 


0.882 




562 


-0.879 ± 0.102 


0.292 ± 0.010 


0.076 


0.267 


0.863 


{MJM,)\M, 


562 


7.635 ± 0.430 


-0.782 ± 0.042 


0.374 


0.474 


-0.652 




562 


9.943 ±0.816 


-0.988 ±0.077 


0.455 


0.450 


-0.515 
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Fig. 1. — Distribution of the g- {green dot-dashed line), r- {orange long-dashed line) and 
i-band {red solid line) SDSS sky levels for the 3041 galaxies in Sample A. 
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Fig. 2. — Comparison of the SDSS and Full sky levels against the manually-computed sky 
measurements for thirty of the largest (black points) and smallest (blue points) galaxies in 
the sample. The labelled points correspond to (a) UGC 7524, (b) UGC 2173, (c) UGC 8334, 
(d) UGC 5882, (e) UGC 7989 and (f) UGC 6346, which show the largest deviations from 
the manually-extracted sky levels in the selected sample. The six galaxies all exceed 7.5' in 
diameter. 
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Fig. 3. — z-band surface brightness profiles for the galaxy UGC 5651 showing the 0.2%, 
0.5% and 1.0% sky flux error envelopes as solid and dotted lines respectively. The fi = 23.5 
mag arcsec"^ surface brightness level is shown as a horizontal dashed line. 




Fig. 4. — Median surface brightness error as a function of surface brightness for all Sample 
A galaxies. The g-, r- and z-band data are shown as green dot-dashed, orange long-dashed 
and solid red lines respectively. 
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Fig. 5. — Sample of nine surface brightness profiles measured by authors MH in blue and 
YZ in red. The solid line shows the exponential disk extrapolation and the disk fit baseline 
is delineated by the green vertical dot-dashed lines. 
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Fig. 6. — Distribution of SDSS-derived axial ratios versus mean g—i galaxy colour (computed 
at the effective radius) for 871 edge-on galaxies as classified by the Galaxy Zoo and ourselves. 
The red dots are median values of the disk flattening, = c/a, in bins oi g — i colours and 
morphological types (upper axis). The histogram on the right includes all 871 galaxies and 
peaks at go = 0.23. Compare with the similar histograms by Lambas et al. (1992; Fig. 6c) 
and Giovanelli et al. (1994; Fig. 7). 
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Fig. 7. — Distribution of physical parameters for 3041 galaxies in Sample A (red) and 652 
galaxies in Sample D (blue). 
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Fig. 8. — SDSS Petrosian radius Rp and apparent magnitude nip j. against the "concentra- 
tion" log(i?p9o— -Rp5o) in the r-band . Red crosses represent galaxies with failed measurements 
of Rp in the SDSS pipeline. The mis-identified targets, which are mostly compact objects, 
lie to the left of the vertical dashed line at log-Rpgo — -Rpso = 0-55. These and red cross 
objects were discarded from our "clean" SDSS sample. 
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Fig. 9. — Comparison of the i-band effective radii and the SDSS effective Petrosian 
radius i?p5o in three inchnation bins: i < 50° {red crosses), 50° < i < 75° {orange triangles) 
and i > 75° {blue circles). Inchnation is the main source of scatter between the Petrosian 
(circular aperture) and isophotal (elhptical aperture) radii. 
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Fig. 10. — Comparison of i-band radial parameters Rp, -Rpso and -Rpgo from our "clean" 
SDSS data sample with our measurements of R23.5, Re- No corrections have been applied to 
the data. The black dashed line has slope unity. 
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Fig. 11. — The comparison of uncorrected i-band radial measurements against the depro- 
jected rotational velocity Kot- Black points show our best Sample D overlaid onto Sample 
A shown in gray. The red dashed line represents the orthogonal linear fit to Sample D. The 
standard deviation a^v and Pearson r correlation displayed in the bottom right corner of 
each panel. The top and bottom panels show Re, R23.5 and R^, and the three Petrosian radii 
for the "clean" SDSS sample, respectively. 
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Fig. 12. — The la standard deviations of all variants of the RV relation plotted against 



rotational velocity in Figure [TT] Measurements are plotted for each parameter in the g, 
r and z-bands in green, orange and red respectively. Sample A is plotted in circles, B as 
triangles, C as squares and D as five-point stars. 
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Fig. 13. — Surface brightness profiles for the galaxy UGC 5651. The isophotal radius, -R23.5, 
corresponds to the location where the surface brightness profile reaches 23.5 mag arcsec"^ 
in each band. These are depicted by the coloured vertical lines. 




Fig. 14. — Comparison of the apparent Petrosian magnitude mp from the "clean" sample of 
SDSS data products and our measurements of mext, 11123,5, and m^ut- The red dashed 
line represents a one-to-one correlation. 
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Fig. 15. — Comparison of the corrected i-band luminosity measurements against the rota- 
tional velocity, Kot- Sample D is shown as black points overlaid on Sample A in gray points. 
The red dashed line represents the orthogonal fit to Sample D. The standard deviation oi^y 
and Pearson r correlation coefficient for the fit to sample D are shown at the bottom of each 
panel. 
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Fig. 16. — The la standard deviations of all luminosity measures compared against rotational 
velocity in Figure 15 Symbols are as in Figure 12 
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Fig. 17. — The 1-a standard deviations of all variations of the RL relation labelled across 
the X-axis with their RL scatters. Measurements in the g, r and i bands are shown in green, 
orange and red respectively with Sample A plotted in triangles, B in circles, C in squares 
and D in five-point stars. 
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Fig. 18. — The VRL relation with Vmt, R23 5 -C/gxt- The red sohd hne is the orthogonal fit 
to the galaxy sub-sample D shown with black points. The red dashed lines are 2-a deviations 
and the gray points show Sample A. The 1-a scatter, Pearson r coefficient and slope of the 
VL {top left), RL {bottom left) and RV {bottom right) relations are shown in the top left 
corner of each panel. The magnitude and direction of the 20% distance uncertainties are 
shown in the lower right corners of each panel. The fitting parameters are listed in Tables 
[9l[9land|9l 
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Fig. 19. — The VRL relation for the galaxy Sample D with distances derived from Vcmb red- 
shifts (black points) and the peculiar velocity corrected distances (red points) from S05/S07, 



with luminosities in the i-band. The line type and parameters are as in Figure 18 The black 



points, fits and line statistics are as in Figure 18 
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Fig. 20. — Dependence of the VRL relation on morphology displayed with both radial scaling 
parameters -R23.5 (left) and (right). The slope and 2-cr deviation of the orthogonal fits 
to Sample D are shown as solid and dotted black lines, respectively. The data points are 
coloured according to Hubble type, as provided by SOS (but not SOT). 
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Fig. 21. — Same as Figure 20 but against colour, {g — i) 
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Fig. 22. — Same as Figure 20 but against surface brightness, fi^. 




8 9 1011 12 8 9 1011 12 



log (M* [Mj) log (M,^^ [Mj) 



Fig. 23. — The stellar mass (top left) and baryonic Tully-Fisher (top right) relations. See 
text about the computation of masses. Ratio of gas mass, Mg, to stellar mass, M*, as a 
function of M* (bottom left) and baryonic mass (bottom right). The slope and 2-a deviation 
of orthogonal fits for Sample D (black points) are shown as solid and dotted red lines, 
respectively. Sample A galaxies are shown in gray. 
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Fig. 24. — Dependence of the stellar mass and baryonic TuUy- Fisher relations on g — i colour. 



